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Abstract 

An  analysis  of  Bolar  velocity  data  obtained  at  the  Stanford  Solar  Observatory 
has  shown  the  existence  of  solar  global  oscillations  (Delache  and  Scherrer,  Nature, 
in  press).  The  oscillations  are  in  the  range  45  to  105  /*Hz  (160  to  370  minutes)  and 
are  interpreted  as  internal  gravity  modes  of  degree  1=1  and  1=2. 
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The  Data 

Observations  of  the  global  solar  velocity  field  have  been  recorded  at  the  Stan¬ 
ford  Solar  Observatory  since  1976.  These  observations  have  recently  been  exam¬ 
ined  for  long  period  oscillations  at  periods  other  than  160  minutes.  Delache  and 
Scherrer  (1963)  have  reported  the  results  of  this  analysis.  This  report  will  review 
that  work  and  expand  on  the  methods  used  to  find  significant  peaks  in  the  power 
spectrum  in  the  presence  of  noise  and  gaps  in  the  data. 

Soherrer  and  Wilcox  (1983)  have  described  the  observing  procedures  and  prel¬ 
iminary  data  reduction  applied  to  the  observations.  The  observations  are 
differential  measures  of  the  llne-of-sight  velocity  of  the  solar  surface.  They  are 
made  by  comparing  the  average  doppler  shift  from  the  center  of  the  solar  disk  with 
the  doppler  shift  from  a  concentric  annulus.  Since  it  appeared  that  interesting 
features  were  present  in  the  low  frequency  part  of  the  spectrum,  the  computation 
was  extended  down  to  i/=0.  Since  only  the  low  frequency  part  of  the  spectrum  is 
examined  here,  the  data  was  averaged  into  5-minute  intervals  and  normalized 
within  each  day  as  in  the  previous  analysis. 

First  the  Fourier  transform  of  the  4  years  of  Stanford  data  (1977-1980)  was  com¬ 
puted  using  a  standard  fast  Fourier  transform  code  (FFT).  It  was  found  that  the 
largest  power  is  in  the  range  45  pH*  -  105  /iHz  (about  160  to  370  minutes).  The 
resulting  power  spectrum  shows  a  number  of  sets  of  lines  with  separations 
corresponding  to  day  side  bands  (i.e.  11.57  /xHz)  which  obviously  come  from  the 
nightly  gaps  in  the  data.  Comparing  the  4-year  spectrum  with  individual  yearly 
spectra,  it  was  found  that  the  strongest  lines  were  more  prominent  in  1979  than  in 
the  other  years.  This  may  be  due  to  the  relatively  clear  skies  and  the  distribution 
of  observing  times  in  1979.  As  a  first  step,  the  analysis  was  restricted  to  that  year. 
Figure  la  shows  the  1979  spectrum.  The  observations  were  begun  on  7  April  and 


continued  through  23  July  with  most  of  the  data  collected  in  late  May  through  July. 
There  were  a  total  of  240  hours  of  data  available.  As  described  in  the  previous 
report  (Scherrer  and  Wilcox,  1063)  the  data  were  seeded  for  each  observing  run  to 
have  the  same  daily  variance.  The  daily  normalization  factor  used  is  highly  corre¬ 
lated  to  the  integrated  signal  power  in  the  5-minute  band,  suggesting  that  the  varia¬ 
tions  in  power  are  due  to  instrumental  changes.  Although  the  analysis  was  done 
using  this  daily  normalized  data,  the  results  were  rescaled  using  the  average  nor¬ 
malizing  factor  to  provide  an  approximate  scale  in  m/s.  This  normalization  has 
both  desired  and  undesired  effects.  It  allows  combining  data  from  different  years 
and  from  different  observation  sets  where  the  raw  data  shows  jumps  in  average  sig¬ 
nal  due  to  variations  in  the  instrument  with  time  and  to  variations  in  the  prelim¬ 
inary  data  reduction  procedures.  This  method  also  removes  apparent  sensitivity 
changes  that  may  result  from  variations  in  solar  activity.  The  main  structure  of  the 
spectrum  is  the  same  for  both  the  normalized  and  raw  data,  so  the  normalized  data 
were  used  to  allow  direct  comparisons  with  different  intervals  of  time. 

Significance  of  Spectrum 

An  estimate  of  the  number  of  peaks  in  the  spectrum  that  are  significantly 
above  the  noise  can  be  made  by  examining  the  cumulative  distribution  of  the  spec¬ 
trum.  A  power  spectrum  computed  from  a  normally  distributed  noise  source  will 
have  an  exponential  distribution,  thus  the  logarithm  of  the  cumulative  distribution 
will  decrease  linearly  with  power.  A  departure  from  a  straight  line  is  an  indication 


Figure  1.  The  power  spectrum  of  velocity  observations  from  the  Stanford  Solar  Observatory 
in  1979.  The  spectrum  in  the  range  45-106  pHs  (300  to  100  minutes)  Is  shown.  Part  (a)  shoos 
the  original  spectrum,  pert  (b)  shews  the  spectrum  to  the  same  sods  after  fourteen  pooka 
wars  identified  and  the  associated  sinusoidal  waves  subtracted  from  the  data.  Tbs  scale 
shown  has  boon  corrected  far  the  average  normalising  factor.  From  Delache  and  Scherrer 
(1963). 
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of  the  presence  of  significant  spectral  features  and  the  slope  of  the  line  is  a  meas¬ 
ure  of  the  variance  (<r)  in  the  spectrum. 

Figure  2  shows  this  plot  for  the  45-105  /iHz  range  of  the  1979  spectrum.  The 
actual  data  is  represented  by  large  dots.  Iha  smaller  points  cure  for  a  spectrum 
computed  the  same  way  from  data  constructed  by  taking  the  original  data  for  each 
day  in  reversed  order,  thus  keeping  the  original  window  function.  Any  coherent  sig¬ 
nal  present  in  the  original  data  will  be  eliminated  by  this  procedure.  The  noise 
statistics  for  the  modified  data  power  spectrum  should  be  unchanged  for  periods  up 
to  the  average  daily  observing  time,  i.e.  tor  frequencies  higher  than  40  /jHz.  It  can 
be  seen  that  in  the  actual  data  there  are  more  peaks  with  value  above  0.25  (m/s)2 
than  in  the  reversed  data.  The  level  at  which  the  actual  distribution  begins  to 
depart  from  the  noise  distribution  is  around  2.5  a.  Since  the  spectrum  was  com¬ 
puted  with  a  resolution  of  0.02  fiHz  but  has  a  natural  resolution  of  only  0.11  fiHz, 
there  will  be  about  five  points  shown  for  each  significant  peak  in  the  spectrum. 
Also,  since  the  data  has  gaps  at  night,  the  day  sidelobe  structure  introduces  2  to  4 
apparently  significant  artificial  peaks  for  each  true  peak.  This  reduces  the  number 
of  independent  spectral  estimates  to  about  200.  These  considerations  suggest  that 
there  are  about  10  independent  peaks  above  the  2.5  a  level. 


Figure  2.  The  cumulative  distribution  function  of  the  speotrum  in  Figure  la. 
The  log  of  the  number  of  spectral  estimates  with  size  larger  than  a  given  value  is 
plotted  vs  that  value.  The  large  dots  correspond  to  the  observed  spectrum.  The 
•mall  dote  are  from  a  similar  spectrum  computed  with  the  data  within  each  day 
taken  in  reversed  order.  The  small  dots  then  refer  to  a  spectrum  with  the  seme 
window  and  noise  characteristics  but  no  coherent  oscillations.  The  variance  which 
Is  deduced  from  the  slope  is  0*s(2Oo m/s)1. 
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Figure  3.  The  spec¬ 
trum  of  the  original 
data  (a),  the  spec¬ 
trum  after  removing 
the  largest  peak  (b), 
after  removing  the 
second  peak  (c),  and 
after  removing  the 
third  peak  (d).  (from 
Scherrer,  1983). 


FREQUENCY  (pHz) 


Table  1.  Frequencies  of  the  14  peaks  found 
in  the  45*105  fi Hz  range  of  the  1079  spec¬ 
trum.  The  frequencies  are  in  /*Hz  with  the 
corresponding  periods  shown  in  minutes. 
The  power  is  shown  as  amplitude  squared 
(m/s)8.  The  classification  of  n  and  1  are 
described  in  the  text.  The  2  peaks  in 
marked  with  an  asterisk  are  an  alternate 
identification  for  the  peak  at  82.20  fiHz  and 
its  associated  side-lobes.  Note  that  the 
peak  at  59.52  /iHz  could  be  identified  as 
part  of  either  an  t-1  or  (=2  series  and  is  in¬ 
cluded  in  the  table  twice. 


CUaning  th e  Spectrum. 

Since  the  spectrum  appears  to  be  dominated  by  the  sidelobe  structure  from 
the  observing  times,  a  procedure  must  be  performed  to  eliminate  these  sidelobes. 
An  iterative  peak  removal  technique  wee  used  to  find  and  remove  the  peaks  one  at  a 
time.  First  an  FFT  was  computed  and  the  largest  peak  determined.  Next  that  peak 
was  accurately  found  with  a  fine  resolution  simple  Fourier  transform  in  the  vicinity 
of  the  peak.  Finally  the  corresponding  sinusoidal  signal  was  subtracted  from  the 
original  data.  This  procedure  was  repeated,  producing  a  list  of  frequencies  free  of 
day  sidelobes.  Figure  3  shows  the  original  spectrum  and  the  spectrum  after  remov¬ 
ing  each  of  the  first  3  peaks.  It  can  be  seen  that  tor  each  peak  removed,  the  entire 
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sidelobe  structure  is  also  removed.  Selecting  peaks  with  amplitudes  greater  than 
about  30  cm/s,  14  peaks  were  found  in  the  frequency  band  45  to  105  fi Hz.  Using  the 
same  amplitude  criterion,  another  8  peaks  were  found  at  frequencies  in  the  range 
extending  to  165  /zHz.  These  higher  frequency  peaks  will  not  be  discussed  here. 
The  list  of  14  peaks  is  given  in  Table  1.  Note  that  since  the  spectrum  is  assumed  to 
consists  of  a  set  of  line  spectra,  the  scaling  used  in  this  paper  is  amplitude  squared 
rather  than  power  in  the  usual  sense.  To  convert  to  the  usual  power  spectrum  scal¬ 
ing  all  values  must  be  divided  by  4  and  would  refer  to  intervals  on  0.1 1  /i Hz  (see  e.g. 
Bath,  1974,  pl83). 

To  see  the  effect  of  removing  the  complete  set  of  peaks  from  the  original  data, 
the  spectrum  of  the  residual  data  after  subtracting  these  22  sinusoids  is  shown  in 
Figure  lb.  Figure  lb  is  shown  to  the  same  scale  as  Figure  la.  It  is  clear  that  far 
more  than  14  peaks  in  the  band  shown  have  been  removed  in  the  resulting  spec¬ 
trum.  This  method  of  peak  identification  appears  to  be  a  powerful  tool  for  analysis 
of  power  spectra  computed  from  data  with  complicated  windows.  The  method  finds 
a  table  of  periods,  phases  and  amplitudes  that  best  clean  the  spectrum.  From 
experiments  with  artificial  data,  we  have  concluded  that  while  the  correct  peak  is 
usually  found,  a  side  lobe  introduced  by  the  window  is  occasionally  selected.  This  is 
a  weakness  of  the  procedure  that  requires  care  in  examining  the  resulting  peaks. 
Ve  have  also  found  that  the  amplitudes  determined  from  data  with  a  complicated 
window  are  not  to  be  relied  upon. 

The  method  described  above  relies  upon  the  single  assumption  that  the  exam¬ 
ined  signals  consist  of  coherent  sine  waves.  It  is  then  not  surprising  that  for  each 
frequency  removed,  a  complete  day-sidelobes  set  of  lines  disappears  from  the  spec¬ 
trum.  A  further  check  is  necessary  to  determine  whether  or  not  the  signal  is 
coherent  tor  the  full  span  of  the  observations.  For  that  purpose,  the  data  with  the 
tabulated  frequencies  removed  was  divided  into  two  parts,  7  April  through  15  June 
and  16  June  through  83  July.  Figure  4  shows  a  set  of  nine  spectra  resulting  from 
this  division.  The  top  row  shows  the  original  spectrum  (left),  the  spectrum  of  the 
first  half  of  the  data  (middle),  and  the  spectrum  of  the  last  half  of  the  data  (right). 
A  large  variation  in  the  appearance  of  the  spectra  is  obvious.  A  first  impression  is 
that  these  spectra  represent  either  very  noisy  data  or  unrelated  data.  Examination 
of  the  bottom  two  rows  shows  that  the  first  impression  is  wrong.  The  second  row 
contains  the  residual  spectra  similar  to  Figure  lb  for  each  of  the  intervals  in  the 
top  row.  The  middle  row,  left  spectrum  is  essentially  the  same  as  Figure  lb.  The 
eenter  and  right  frames  are  for  two  halves  of  the  interval  but  cleaned  by  subtract¬ 
ing  the  sinusoids  found  for  the  entire  interval.  The  spectra  of  each  part  separately 
da  not  show  any  of  the  removed  frequencies.  If  any  of  the  frequencies  removed 
were  net  present  in  both  halves  of  the  data,  they  would  have  been  artificially  put 
into  the  cleaned  data  by  the  subtraction  of  the  sinusoid  over  the  entire  interval.  In 
that  case  they  would  show  in  the  spectrum  of  one  or  both  halves  of  the  cleaned 
data.  This  is  strong  evidence  that  the  signals  are  all  present  in  both  halves  of  the 
observations.  The  bottom  row  demonstrates  the  variability  in  amplitude  caused  by 
the  observing  window.  Using  only  the  times  of  observations  and  the  list  of  frequen¬ 
cies  and  phases  found  above,  a  cleaned  dataset  was  reconstructed.  The  spectra 
from  the  entire  interval  and  the  two  halves  are  shown  in  the  bottom  row.  Note  that 
all  throe  frames  are  computed  from  artificial  data  generated  from  the  same  set  of 
frequencies  and  phases.  The  only  difference  between  frames  on  the  bottom  row  is 
the  observing  window.  The  conclusion  is  that  for  windows  obtained  at  Stanford 
(which  are  similar  to  the  average  distribution  of  observations  in  the  Crimea  and  at 
SCLERA),  the  relative  amplitudes  of  peaks  may  be  dominated  by  interference 
caused  by  data  gaps. 
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Peak  Identification 

Returning  to  table  1,  two  of  the  14  peaks  are  apparently  day  aliases  at  1/4  and 
1/B  of  a  day.  (With  only  3  months  of  data  the  previously  reported  peak  at  160.01 
minutes  can  not  be  distinguished  from  1  /B  of  a  day  and  is  masked  by  the  sidelobes 
of  the  lower  day-harmonics).  This  leaves  12  peaks,  for  further  analysis.  There  is 
some  ambiguity  in  the  true  identification  of  the  largest  of  the  remaining  peaks.  The 
first  1/day  sidelobe  is  of  the  same  amplitude  as  the  main  peak.  This  could  be  due 
to  two  true  peaks  with  separation  1/day  or  2/day.  By  careful  examination  of  the 
peak  shapes,  the  peak  at  62.29  /xHz  has  been  identified  as  the  true  peak  with  a 
smaller  peak  at  96.94  fMx,  although  the  identification  could  have  been  made  as  a 
true  peak  at  73.64  /xHz  with  smaller  peaks  at  50.72  and  96.94  /*Hz.  The  alternate 
identifications  are  shown  in  parentheses  in  Table  1. 

Interpretation 

In  order  to  interpret  the  peaks  found  in  this  spectrum,  one  must  seek  guidance 
from  theory.  Estimates  for  the  spacing  of  g-mode  oscillations  made  from  the  full 
asymptotic  approximation  (Tassoul  I960)  and  calculated  from  complete  solar 
models  (Berthomieu,  Provost,  and  Christensen-Dalsgaard,  private  communications) 
suggest  that  the  g-mode  oscillations  of  the  same  degree  l  should  be  about  equally 
spaced  in  period  for  order  large  enough  (n  >  6).  Therefore  the  list  of  prominent 
peaks  was  examined  for  equal  spacing  in  period.  Three  of  the  four  largest  peaks  are 
separated  in  period  by  15.5  minutes.  This  is  very  close  to  the  spacing  of  the  stan¬ 
dard  solar  model  for  g-modes  with  degree  two  and  order  greater  than  10. 

To  aid  in  further  identification  of  the  modes  the  asymptotic  approximation 
(Tassoul  I960)  is  used.  This  approximation  shows  that  for  sufficiently  large  order  n 


figure  5.  Tbs  observed  modes  are 
plotted  within  the  structure  of  tbs 
asymptotic  formula  described  in  tbs 
text.  The  line  Is  determined  from 
tbs  three  largest  peeks  that  have 
been  assigned.  Prom  Delache  and 
Se  horror  (1963). 


7 


the  period  T  is: 

T  =  Tm  («  +  g~  ^  /  vTTTTT) 

If  the  three  peaks  with  15.5  minute  separation  are  part  of  an  1=3  series,  then  T9  and 
thus  the  probable  spacing  for  the  1=1  and  1=3  series  can  be  computed.  Doing  this, 
it  is  seen  that  most  of  the  remaining  peaks  are  likely  to  be  part  of  the  1=1  series. 
Using  the  estimate  for  1,  the  observed  periods  can  be  plotted  as 
T  Vt(l  +  l)  vs  (k  +  1/ 2  -  1/  4)  where  k  is  an  integer  increasing  with  period.  If  the 
peaks  are  consistent  with  the  model  and  the  correct  degree  1  has  been  assigned, 
the  points  will  all  lie  on  a  straight  line.  From  the  intercept  of  this  line  both  the  best 
value  for  n  and  the  correct  value  for  -1/4  term  can  be  determined.  This  has  been 
done  in  Figure  5. 

The  three  largest  peaks  were  used  to  define  the  1=2  series  and  thus  to  find  T0 . 
The  rightmost  dots  in  Figure  5  represent  the  three  largest  peaks  which  determine 
the  1-2  series.  The  other  modes  identified  in  Table  1  are  also  shown.  The  line  was 
found  from  the  three  large  1=2  peaks  only.  It  can  be  seen  that  the  other  peaks  are 
organized  in  the  expected  g-mode  structure.  The  T9  implied  by  the  identified 
modes  is  38.0*0.5  minutes.  The  intercept  of  the  best  fit  line  with  the  abscissa  is 
expected  to  be  0.35  and  is  found  to  be  0.2*0. 1,  thus  the  order  n  is  most  likely 
correctly  determined.  If  n  is  decreased  by  one  for  each  peak,  T9  becomes  about  41 
minutes  but  the  scatter  is  increased  somewhat. 

The  largest  peak  at  63.39  /zHz  appears  not  to  be  consistent  with  the  asymptotic 
formula  with  the  T9  we  have  found.  If  the  alternate  choice  of  the  63.39  p.Hz-73.84 
/zHz  pair  were  chosen,  the  largest  peak  would  be  at  73.84  /zHz  which  would  be 
identified  as  1=1,  n=B.  In  this  case  a  peak  at  50.7  /zHz  shows  up  in  the  peak  finding 
procedure  and  would  be  identified  as  1=2,  n=30. 

Note  that  in  several  cases  two  peaks  have  been  assigned  to  one  order  n  in  the 
1=1  series.  These  peaks  have  an  average  separation  of  1.3  /zHz  and  could  be  evi¬ 
dence  of  rotationally  split  modes.  Assuming  the  1=1  modes  are  m=±l,  the  1.3/tHz 
separation  would  correspond  to  a  rotation  rate  n=7.5xl0-asec_1.  The  rotational 
splitting  kernel  for  these  modes  has  a  broad  maximum  between  0.3  and  0.5  R. 
Assuming  a  smooth  form  for  0(r)  as  in  Figure  3  of  Gough  (1983),  J z  can  be  calcu¬ 
lated  and  is  of  order  1.7xl0“a.  This  value  differs  from  that  found  by  Hill  et  al. 
(1983). 

The  sensitivity  of  the  Stanford  instrument  to  degree  1=2  g-mode  oscillations 
must  be  considered.  For  the  previously  reported  acoustic  mode  oscillations  in  the 
5-minute  range,  the  Stanford  velocity  differencing  scheme  is  most  sensitive  to 
modes  of  degree  1=3  to  1=5  but  is  not  sensitive  to  modes  of  degree  1=1  or  1=3.  For 
low  frequency  g-mode  a,  Gough  and  Cristensen-Dalsgaard  (1983)  have  shown  that  for 
periods  around  160  minutes  the  Stanford  instrument  is  most  sensitive  to  modes 
with  degree  1=5  to  1=7.  For  longer  periods  however,  the  instrument  begins  to  be 
relatively  more  sensitive  to  modes  with  lower  degree.  This  variation  with  period, 
unlike  the  case  for  the  p-modes,  is  due  to  the  importance  of  the  horizontal  motions 
induoed  by  the  internal  gravity  modes.  To  properly  estimate  the  sensitivity  to  the 
various  modes  the  oscillation  energy  must  also  be  considered.  In  the  simplest  case, 
the  combined  effect  of  equlpartition  of  energy  between  modes  and  the  mode  visibil¬ 
ity  with  the  Stanford  instrument  yields  the  greatest  sensitivity  for  1=1  and  1=3  at 
these  longer  periods. 
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Figure  6.  Artificial  spectra  computed  from  a  single  frequency  signal  with  a  selec* 
tion  of  sample  times  and  lifetimes.  The  sample  duration  T  and  mode  lifetime  r  vary 
from  3  days  to  100  days. 
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The  theory  predicts  several  hundred  periods  (or  low  degree  g-modes  in  the 
period  interval  examined.  Only  the  most  visible  handful  of  peaks  has  been  exam¬ 
ined  and  has  been  found  to  be  consistent  with  the  theory.  The  relative  peak  sizes 
and  simplicity  of  the  analysis  leads  to  the  present  mode  identification.  Alternative 
identifications  are  quite  possible  but  they  can  only  be  tested  with  more  data  from 
other  years  or  from  other  observatories. 

Note  Concerning  Lifetimes 

There  are  some  general  considerations  concerning  the  appearance  of  power 
spectra  of  data  containing  oscillations  of  finite  lifetime  that  must  be  remembered 
when  interpreting  such  spectra.  It  is  generally  recognized  that  the  spectrum  of  an 
oscillator  with  a  lifetime  r  greater  than  the  observation  time  T  will  be  a  "line"  that 
will  appear  as  a  vine8  with  width  1/  T.  It  is  also  recognized  that  an  observation  of  a 
single  oscillation  event  with  lifetime  t  and  observation  time  T  will  result  in  a  similar 
spectral  line  but  with  width  1/  r.  What  is  sometimes  forgotten  is  that  if  the  observa¬ 
tion  duration  includes  severed,  but  not  memy  instances  of  an  oscillation  lifetime,  i.e. 
T> t  but  not  T»t.  the  spectrum  will  be  a  complex  structure  of  lines  of  width  1/  T 
in  an  envelope  of  width  1/t.  If  T  ~2r  the  spectrum  will  look  similar  to  that  of  a 
rotationally  split  1  =  1  mode.  Figure  6  shows  spectra  computed  from  artificial  data 
constructed  with  T  =  3-days  through  T  =  100-days  and  r  =  3  days  through  t  »  100 
days.  There  were  no  data  gaps  assumed.  The  plot  scales  were  adjusted  such  that 
the  plots  almost  fill  the  frame.  An  examination  of  Figure  8  shows  that  caution  must 
be  used  in  interpreting  spectra  of  oscillatory  signals  of  unknown  lifetimes. 
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